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ABSTRACT 

We present a study of available Hubble Space Telescope (HST) spectroscopic and 
imaging observations of the low ionization nuclear emission line region (LINER) galaxy 
NGC 1052. The WFPC2 imagery clearly differentiates extended nebular Ha emission 
from that of the compact core. Faint Object Spectrograph (FOS) observations provide 
a full set of optical and UV data (1200-6800 A). These spectral data sample the in- 
nermost region (0."86 x 0."86 ~ 82pc x 82pc) and exclude the extended Ha emission 
seen in the WFPC2 image. The derived emission line fluxes allow a detailed analysis 
of the physical conditions within the nucleus. The measured flux ratio for Ha/H/3, 
Fffoi/Ffl^ = 4.53, indicates substantial intrinsic reddening, E(B-V)=0.42, for the nu- 
clear nebular emission. This is the first finding of a large extinction of the nuclear 
emission line fluxes in NGC 1052. If the central ionizing continuum is assumed to be 
attenuated by a comparable amount, then the emission line fluxes can be reproduced 
well by a simple photoionization model using a central power law continuum source 
with a spectral index of a = —1.2 as deduced from the observed flux distribution. A 
multi-density, dusty gas gives the best fit to the observed emission line spectrum. Our 
calculations show that the small contribution from a highly ionized gas observed in 
NGC 1052 can also be reproduced solely by photoionization modeling. The high gas 
covering factor determined from our model is consistent with the assumption that our 
line of sight to the central engine is obscured. 

Subject headings: galaxies: active — galaxies: individual (NGC 1052) 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data archive at the Space 
Telescope Science Institute. STScI is operated by the Association of Universities for Research in Astronomy, Inc. under 
NASA contract NAS 5-26555. 



1. Introduction 

Low ionization nuclear emission line regions, 
(LINERs) have been found to occur in nearly 1/3 
of pll bright, nearby galaxies j Ho eL al. 1995 ): — 
They share many of the emission characteris- 
tics of active galactic nuclei (AGN), but are dis- 
tinguished by their low luminosities and rela- 
tively strong, low- ionization emission lines (e.g., 
[O I] AA6300, 6363, [O II] A3727, [N II] AA6548,6584, 
[S II] AA6716,6731). A recent study by Ho (1999) 
has demonstrated that the UV — X-ray spec- 
tral energy distribution of LINERs is harder than 
in higher luminosity AGNs. This is consistent 
with the general trend exhibited in AGNs of de- 
creasing a ox with decreasing luminosity (Avni 



nated by low ionization features. For example, 
the [O I] A6300 and [O II] A3727 fluxes exceed 
that of [O III] A5007, plus [N II] A6584/ Ha >1 
([Fosbury et al. 1978; , [Ho et al. 1993j , |Ho et al 



Tananbaum 1982). The physical relation be- 
tween LINERs and AGNs is unclear. They may 
be simply the low luminosity extension of the 
AGN phenomenon, powered by accretion onto a 
compact object or, alternatively, they may be the 
result of entirely different mechanisms. 

Although central source photoionization ap- 
pears to be the dominant mechanism powering 
the nebular emission seen in the high luminos- 
ity Seyfert galaxies, it has not been thoroughly 



tes 



ed fui galaxies cunlaiiiiiig LINERs. Cumpel- 

ing ionization mechanisms that have been used to 
successfully explain at least some of the emission 
characteristics of LINERs include photoioniza- 
tion by hot stars ([Terlevich &: Melnick 1985 , Fil- 



ippenko fc Terlevich 1985| ), photoionization by 



non-thermal emission typical of luminous AGN 



(Ferland & Netzer 1983|, 


Ho et al. 19931) and 


shock-heating (Heckman 1980, 


Dopita & Suther- 


land 1995, 


Dopita 1996|). 


LINERs might repre- 



sent a heterogeneous class of objects since the 
proposed ionization mechanisms can all explain 
the low apparent luminosity of LINERs. Perhaps 
more than one of these various mechanisms is at 
work in the same object. 

NGC 1052 is a nearby (z=0.0045; |de Vau 



ouleurs et al. 1991) elliptical galaxy (E4) with 



1997|) . Fosbury et al. 1978, hereafter F78, showed 
that the emission line gas is comprised of a lu- 
minous central component (d < 2") and a dif- 
fuse extended region (d~20"). Each emission 
region contributes approximately equally in to- 
tal Balmer line flux. Furthermore, the extended 
component emits many of the other prominent 
optical lines that are observed in the spectrum 
of the central region, such as [O II] A3727 and 
[O III] A5007 (F78). 

The observed optical continuum emission is 
dominated by the host galaxy, with no clear de- 
tection of any underlying non-thermal emission 



(F78; Fosbury et al. 1981 ). Radio studies, on the 
other hand, clearly show it to have both compact 
and extended continuum emission components. 
At ~1.5 GHz, its radio emission is characterized 
by a flat, unresolved (< 1") core and two weaker 
lobes extending 10-15" on either side of the com- 
pact source QWrobel 1981 ). NGC 1052 has been 
shown to have an IR excess in the 5-20 /jm band 
within the inner 2" of the nucleus flBecklin~et 



al. 1982| ). A recent analysis of ROSAT/RRI im- 
agery shows some extended soft X-ray emission, 
but the low spatial resolution at these energies 
prevents isolating any compact emission at scales 
less than 10" ( |Guainazzi fc Antonelli 1999| ). X- 
ray observations with ASCA indicate a very flat 
spectrum for the hard X-ray emission (2-10 keV) 
of NGC 1052 (|Guainazzi fc Antonelli 19991 ). 



a strong nuclear emission line spectrum domi- 



As with LINERs in general, there is much 
disagreement on the nature of the excitation 
mechanism responsible for the nuclear activity 
in NGC 1052. Initial findings have suggested 
that shock-heating could produce the observed 
line emission QKoski fc Osterbrock 1976| ; F78). 
This claim is largely based upon the high elec- 
tron temperature (T~33,000 K) implied by the 
strong [O III] A4363 emission. Furthermore, sim- 
ple shock excitation models have been shown to 
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match the prominent optical emission features 
( Koski fc Ostcrbrock 197*61 ; F78). In an alternate 
interpretation, Ferland & Netzer (1983), and 
more recently Ho et al. (1993), have shown that 
many of the strong emission lines observed in 
NGC 1052 (and in LINERs in general) could be 
produced by photoionization from a nonthermal 
continuum typical of Seyfert galaxies, namely a 
power law with spectral index, a~-1.5. The ma- 
jor difference in the models between the typi- 
cally more luminous Seyferts and the LINERs 
is that the ionizing flux is inferred to be sig- 
nificantly weaker in LINERs. Subsequent de- 
tailed photoionization modeling of NGC 1052 
by Pequinot (1984), hereafter P84, has shown 
that many of the optical emission lines, includ- 
ing [O III] A4363, could be fit well by photoion- 
ization of a multi-density component gas. His 
study assumed a central photoionizing source 
consisting of a black body continuum with an X- 
ray tail extending to higher energies. Recently, 
a polarimetric study has shown NGC 1052 to 
have a broad Ha emission line component with 



FWHM~5000 km s" 1 ( [Barth et alTT^ ). This 



provides compelling evidence for the presence of 
a compact object and lends strong support to the 
non-stellar photoionization hypothesis. 

In this paper, we present a detailed analy- 
sis of the nuclear activity in the prototypical 
LINER galaxy NGC 1052. Specifically, we use 
Hubble Space Telescope (HST) observations of 
NGC 1052 to explore if photoionization repre- 
sents a viable option in explaining the emission 
line spectrum in this LINER. In Section 2 we de- 
scribe the HST imagery and spectroscopic obser- 
vations used in this study and present our mea- 
sured emission line fluxes and an analysis of the 
WFPC2 Ha imagery. We examine the possible 
role of extinction on the emission line fluxes and 
ionizing continuum. The methodology and re- 
sults of our photoionization modeling are given 
in Section 3. In Section 4, we discuss the im- 
plications of our results and compare them with 
previous studies to form a comprehensive picture 



of the activity of NGC 1052. Finally, we summa- 
rize our results and present our conclusions in 
Section 5. 

2. Observations and Emission Line Mea- 
surements 

2.1. FOS Observations 

The spectral data used in this study were 
obtained with the Faint Object Spectrograph 
aboard the Hubble Space Telescope (HST/ FOS) 
on 1997 January 13. They were retrieved from 
the Space Telescope Science Institute (STScI) 
data archives. All observations were of the ac- 
tive nucleus of NGC 1052 through the 0".86-pair 
square aperture (post-COSTAR), covering 1200- 
6800 A. Thus, the aperture sampled the central 
82 pc x 82 pc of the nuclear region of NGC 1052, 
assuming D = 19.6 Mpc for H D = 75 km s _1 Mpc -1 . 
The optical and near-UV data were obtained 
with high resolution gratings giving A/AA~1300, 
while the UV spectrum has a lower resolution 
given by A/AA~250. The individual spectra are 
described in Table 1 along with their exposure 
times, grating settings, and spectral resolutions. 
Because of a guide-star lock failure, the data may 
have slightly degraded spectral resolution and 
sample a slightly larger spatial region than the 
aperture size would indicate. If the extended re- 
gions emit differently than the nucleus, this may 
have an effect on the observed emission line ratios 
and absolute fluxes. 

Using the maps of Burstein &: Heiles (1982), 
we find a foreground Galactic reddening value of 
E(B-V) = 0.02 in the direction of NGC 1052. 
Figure 1 shows the spectrum in the rest frame of 
NGC 1052, corrected for Galactic extinction us- 
ing the extinction law of Savage & Mathis (1979). 

2.1.1. UV Emission 

Examination of the UV spectrum below 2000 A 
in Figure 1 shows a weak UV continuum com- 
pared to the optical continuum, although, there 
appears to be a slight rise toward lower wave- 
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lengths. This UV continuum shows no evidence 
for recent starburst activity in the nucleus of 
NGC 1052. Specifically, the common signatures 
of hot O and B stars, namely broad wind features 
from C IV A1549, Si IV A1400, and N V A1240 
and narrower photospheric lines (Neubig & Bruh- 
weiler 1997, 1999), are absent within the low de- 
tection limits of the data. Although the results 
are not definitive, the absence of stellar features 
is consistent with the UV continuum being pro- 
duced by a non-thermal compact source. Numer- 
ous UV emission lines are detected in the FOS 
spectrum, compared to the IUE study in which 
only two lines, C III] A1909 and C II] A2326, were 



unambiguously measured (Fosbury et al. 1981) 



2.1.2. Emission-Line Measurements 

Visual inspection of the optical and near-UV 
data at A > 2000 A in Figure 1 clearly shows that 
the host elliptical galaxy dominates the contin- 
uum. Strong stellar absorption features are ev- 
ident and may affect flux measurements of cer- 
tain emission lines, especially the Balmer lines. 
In addition, small wavelength variations in the 
composite stellar spectrum of the host galaxy 
may mimic emission line features. To remove the 
spectral contribution of the host galaxy, we have 
used template spectra of elliptical galaxies from 
the HST/FOS archival database to approximate 
its integrated stellar spectrum. We selected a 
proper template spectrum for each of the three 
optical gratings (G270H, G400H, and G570H) 
based upon the following criteria: 
1) We required the prominent absorption fea- 
tures that do not coincide with emission fea- 
tures in NGC 1052 to match those in the host 
spectrum of NGC 1052 (namely, Mg I A2852, 
Ca II K A3934, G Band A4301, Fe I A4440, 
Fe I A4532, Mg lb A5176, Fe I A5333, Ca 1/ 
Mg H A5590, Na I Doublet A5893, and Ca 1/ Fe I 



A6161; Pritchet fe van den Bergh 1976 , Crenshaw 
^ Peterson 1985| ); 

2) We required the overall galactic flux distribu- 
tion to be similar to that in NGC 1052. 



The template spectra were then scaled in mag- 
nitude, so that the strengths of the absorption 
features matched that of NGC 1052, and sub- 
tracted from it. Our template spectra are listed 
in Table 1 and plotted with the corresponding 
spectra of NGC 1052 in Figure 1. Since there is 
very little stellar contribution to the UV in el- 
liptical galaxies and no evidence for hot stars in 
NGC 1052, no template was used for the G160L 
grating. 

We then measured the emission line fluxes us- 
ing the subtracted spectrum corrected for Galac- 
tic extinction. Flux ratios relative to H/3 are 
given in Table 2. To measure the line fluxes, 
uncontaminated continuum regions were selected 
on either side of the features to be measured. A 
linear fit was made to the continuum underly- 
ing the emission features and the lines were ex- 
tracted. In several cases, it was difficult to dis- 
tinguish between real emission lines and artifacts 
remaining from our attempt to remove the under- 
lying spectum of the host galaxy. This is most 
notable in the He II A4686 emission. In Table 2, 
we have indicated where these uncertainties are 
large. The emission line widths are broader than 
typically observed in the narrow line regions of 
AGN, corresponding to a range of velocities be- 
tween 900 km s- 1 < FWHM < 1100 km s" 1 (by 
comparison, the instrumental resolution for the 
high dispersion gratings is FWHM~230 km s _1 ). 

In several cases, it was necessary to deblend 
line contributors before accurate emission line 
fluxes could be derived. In all cases, the spectra 
were converted to velocity space to preserve the 
correct line widths. For example, Ha is heavily 
blended with the [N II] AA6548,6584 doublet. To 
extract the Ha feature, the normalized line pro- 
file of H/3 was used to make a template for both 
of the [N II] lines, since no broad wings were de- 
tected on the H/3 profile. This technique proved 
more accurate than fitting with Gaussian profiles 
because it approximates the actual structure of 
the emission lines. The [N II] template was cre- 
ated by separating the two profiles according to 
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their intrinsic transition wavelengths and assum- 
ing the flux ratio corresponds to that of their Ein- 
stein A-values (A,x6584/-A-A6548 = 2-96). The wave- 
length offset and relative flux of the template 
were then varied until the residual of the [N II] 
contribution subtracted from the Ha±[N II] line 
blend gave the best fit to Ha, namely, a single 
emission feature on a flat continuum (Figure 2a) . 
Figure 2b shows the resulting Ha fit compared to 
the H/3 velocity profile and demonstrates the re- 
liability of this technique. 

Several other line blends were treated using 
the same technique. It is important to de- 
blend and provide accurate flux ratios for the 
[S II] AA6716,6731 lines because they provide a 
low density constraint for the emission line re- 
gion. To fit them, H/3 profiles were centered 
at the intrinsic wavelengths of the two [S II] 
lines. The flux ratio was varied until residu- 
als were minimized. For line blends of other 
hydrogen lines, namely H7+[0 III] A4363 and 
He+[Ne III] A3967, the H/3 profile was used 
template for the Balmer line since it will have an 
identical intrinsic shape. All lines which required 
deblending are noted in Table 2. 

2.1.3. Uncertainties 

The primary sources of uncertainty are con- 
tinuum placement in extracting emission lines 
and the approximations associated with the de- 
blending procedure. For the first case, the un- 
certainty depends on the strength of the emis- 
sion line relative to the noise of the data and 
the existence of host galaxy features. We exper- 
imented with continuum levels for lines of differ- 
ent strengths. For the optical and near-UV data, 
where the host galaxy is dominant but spectral 
resolution is high, this leads to an estimate of the 
uncertainty for "strong" lines of ±5% (defined as 
having F > Fup), for "moderate strength" lines 
±10% (0.2 xF^ < F < F H p), and for "weak" 
lines ±20% (F < 0.2xF H/3 ). For wavelengths 
below ~2000 A, there is much less contribu- 
tion from the host galaxy, however, the signal- 



to-noise ratio is significantly deteriorated. Con- 
sequently, we find flux uncertainties of the mod- 
erate strength lines to be ±20% and the weak 
lines to be ±30%, for wavelengths below 2000 A. 

Additional uncertainties are introduced for 
lines that require deblending. To reliably esti- 
mate this uncertainty, the flux of the residual 
associated with the [N II] A6584 line (see Fig- 
ure 2a) was measured and shown to be ±10% of 
the measured flux for the line. This represents an 
upper limit to the general deblending uncertainty 
because it was the largest relative residual of all 
procedures. These two uncertainties were then 
added in quadrature to obtain the total error for 
the measured line ratios. The total uncertainty 
estimate for each emission line is listed in Table 2. 

2.1.4- Reddening 

To estimate the internal reddening of NGC 
1052, we compared our measured Ha/H/3 flux ra- 
tio (Fffa/^Hp = 4.53, corrected for Galactic ex- 
tinction) to an assumed intrinsic ratio (Ha/H/3)/. 
For N e = 10 4 cm" 3 and T e = 10,000 K and as- 
suming case B recombination, (Ha/H/3)/ = 2.86, 
which is consistent with the typical ratio de- 
termined for AGN (Osterbrock 1989). Using 
the extinction curve of Savage & Mathis (1979), 
this implies a reddening of E(B-V) = 0.42 ±.20, 
where the uncertainty was derived from the er- 
rors associated with the Ha and H/3 flux mea- 
surements. The deduced intrinsic line ratios are 
listed in Table 2. Uncertainties were estimated 
from propagation of errors associated with the 
reddening and observed flux ratios. Since all pre- 
vious studies of NGC 1052 determined negligible 
internal reddening, it is important to examine 
this issue in more detail (Koski & Osterbrock 
1976; F78; Fosbury et al. 1981; Ho et al. 1993; 
Ho et al. 1997). 

The intrinsic He II A1640/He II A4686 ratio is 
theoretically a better measure of reddening than 
Ha/H/3 because of the large span in wavelength 
between the two lines ( |Kraemer et al. 1994] ) . 
However, the large flux uncertainties associated 
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with the He II lines greatly diminish the useful- 
ness of this diagnostic for our data set. Yet, we 
can still use it as a check on our reddening de- 
duced from the Balmer line ratio. We assume an 
intrinsic ratio of (He II A1640/He II A4686) / = 6.6, 
corresponding to N e = 10 4 cm~ 3 and T e = 10,000° K 



flSeaton 1978 ). This is consistent with the dered- 
dened ratio of He II A1640/He II A4686 > 4.2 con- 
sidering the large uncertainties associated with 
these line fluxes (Table 2). 

We can also use the relative flux ratios of the 
higher order hydrogen Balmer lines beyond H/3 as 
a consistency check. A comparison of the dered- 
dened Balmer ratios with their intrinsic values 
is given m Table 3 flOsterbrock 198SQ . They arc 
seen to be consistent within the quoted uncer- 
tainties. 

The colors of the host galaxy of NGC 1052, 
B-V = 0.91, U-B = 0.42, are in accord with typ- 
ical colors found in normal elliptical galaxies of 
similar Hubble type, E4 (F78). Thus, the large 
extinction found above is not consistent if ap- 
plied to the entire galaxy. If the implied red- 
dening is correct, it suggests that dust is local- 
ized to the nuclear region of NGC 1052. There 
is evidence which supports this. The B-R and 
B-I images of Carter et al. (1983) and B-R im- 
age of Davies et al. (1986) clearly show a diffuse 
dust lane extending out from the central region to 
about 20". Furthermore, the IR measurements of 
Becklin et al. (1982) indicate a large IR excess be- 
tween 5-20 fan confined within 2" of the nucleus. 
The strength of the IR flux at these wavelengths 
is not consistent with stellar radiation emanat- 
ing from cool stars and may be thermal radiation 
from heated dust. 

2.2. Ha Imagery 

Images of the nuclear region of NGC 1052 were 
taken with WFPC2/#ST using the F658N fil- 
ter on 1996 June 16. These data were obtained 
from the STScI data archives. This filter iso- 
lates any extended Ha emission. In Figure 4, we 
have combined two exposures (with 800 and 900 s 



exposure times) and removed cosmic ray events. 
The size of the aperture used for the FOS ob- 
servations is shown over the bright emission core 
for reference, although the precise position of the 
aperture is uncertain due to the guide star lock 
failure. Strong filamentary nebular Ha emission 
is seen to extend about 1" around the compact 
nucleus, with a more diffuse halo extending to 
further distances. At a position angle of ~235° 
is a narrow filament of Ha emission. This is com- 
pared to the radio jet/lobe observed by Wrobel 
et al. (1984) at 1.5 GHz which is centered at a 
position angle of ~275°. 

3. Photoionization Models 

The purpose of the following analysis is to de- 
termine if the emission in the compact nuclear 
region of NGC 1052 is consistent with photoion- 
ization by a simple non-thermal continuum, char- 
acteristic of AGNs. We seek the simplest model, 
making the fewest number of assumptions, that 
gives a satisfactory fit to the dereddened intrin- 
sic emission line spectrum. Discrepancies be- 
tween our model results and the observed spec- 
trum then provide further insight into the actual 
physical conditions in NGC 1052. 

Photoionization models were computed using 



Version 90 of the computer code CLOUDY (|Fer 



land 1996| ) . These models were parameterized by 
the gas density n#, continuum shape, and ioniza- 
tion parameter 

Q 



u 



where Q is the number of hydrogen ionizing pho- 
tons s _1 and R the distance from the inner face of 
the cloud to the central continuum source. When 
a good fit was achieved, observational constraints 
on the ionizing luminosity provided constraints 
on R. 

We adopted the following idealized properties 
of the emission line gas. In cases where more 
than one density component was required, we as- 
sumed that each component was photoionized by 
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unattenuated radiation from the central source. 
Thus we assume, a priori, relatively small cover- 
ing factors for any inner gas components. Phys- 
ically, this corresponds to a clumpy gas com- 
posed of radiation bounded clouds. We assume 
that the observed radiation is emitted from the 
"front" faces of the clouds, that is, the side of 
the clouds facing the central continuum source. 
Each component was assumed to have constant 
density. We began by omitting dust grains from 
the emission line gas, hence, in our initial model, 
we assume solar abundances and that all red- 



dening occurs external to the line emitting re- 



gion. Thus, the abundances relative to hydro- 
gen are: He = l.OOxlO -1 , C = 3.55xl0" 4 , 
N = 9.33xl0" 5 , = 7.41xl0~ 4 , Ne = 1.17 x 10" 4 , 
Mg = 3.80xl0" 5 , Si = 3.55X10" 5 , S = 1.62xl0~ 5 , 
and Fe = 3.24xl(T 5 QGrevesse fe Anders 19891 ). 
The effects of grains mixed with the gas will be 
treated in a later section. 

3.1. Choice of Model Parameters: 

We adopted a simple power law ionizing con- 
tinuum parameterized by spectral index a. Ob- 
servations in the UV and X-ray were used to con- 
strain its shape (Figure 3). We assumed that the 
observed UV continuum flux was reddened by the 
same amount as the emission line fluxes, E(B- 
V) = 0.42, as deduced from the Ha/H/3 ratio. 
This gives a projected intrinsic flux at the Lyman 
edge of Fiifyd~1.9xl0 -27 ergs s _1 cm -2 Hz -1 . 
Results from ASCA and ROSATshow NGC 1052 
to have a rather flat spectrum between 2—6 keV, 
a steeper continuum between both 6—10 keV and 
~1— 2 keV, and a flattening at energies below 
1 keV (0 

uamazzi fe Antonelli 1999| ). This indi- 
cates that the soft X-ray emission is likely heavily 
absorbed, which is consistent with the modeling 
by Weaver et al. (1998). If this absorption is 
due to a dusty gas, then our previous assump- 
tion that the UV continuum is attenuated by a 
column of dust is justified. We chose the ob- 



intrinsic nuclear X-ray flux, as shown in Fig- 
ure 3. Interpolating between the estimated in- 
trinsic UV continuum flux at 1 Ryd and the X- 
ray flux at 6 keV gives a = —1.2 and a total 
intrinsic luminosity for the ionizing continuum 
of L>i^ y( i~2.0xl0 42 ergs s -1 . We note the large 
uncertainty in this value due to the uncertainty 
in the amount of obscuration by dust. Previ- 
ous studies have indicated that the continuum 
emission in AGNs is typically not reddened as 
much as the emission line flux, thus, we may be 



overestimating the total luminosity (Costero & 
Qsterbrock 19771). 



In Figure 3 we have plotted the spectral en- 
ergy distribution of the continuum emission from 
the radio through hard X-ray energies of the cen- 
tral point source. In each energy band, we have 
plotted only the emission that cannot be spa- 
tially resolved. Data for both the observed and 
reddening corrected UV continuum are included, 
illustrating the large effect that our assumed ex- 
tinction has on the energetics of the active nu- 
cleus. 

Simple model calculations using a gas with 
constant density do not match the intrinsic emis- 
sion line spectrum of NGC 1052 for any choice of 
density. One important problem with these mod- 
els is that they cannot match both the auroral 
and nebular line fluxes simultaneously. When a 
low density gas is used, n^~10 3 ' 5 cm -3 , the neb- 
ular lines are fit well (i.e., [O III] AA5007,4959, 
[N II] AA6548,6584, [O II] A3727, and [S II] 
AA6716,6731), but their auroral line counter- 
parts are under estimated, namely [O III] A4363, 
[O II] A2470, and [S II] AA4069,4076. For higher 
densities, the opposite is true. As noted ear- 
lier, it was the high temperature implied by the 
[O III] A4363/[0 III] AA4959,5007 ratio that led 
to the original interpretation of shock excitation 
for NGC 1052. Indeed, our [O III] line flux ratio, 
[O III] A4363/ [O III] AA4959,5007~ 0.089, un- 



served ASCA flux at 6 keV, F, 



6keV 



1.4x10- 



-"30" 



plies T>100,000 K in the low density limit (Os 
terbrock 1989), which is a factor ~5 higher than 



ergs s 1 cm 2 Hz , as the best estimate of the 



we would expect in the ++ zone of a photoion- 
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ized nebulae. However, choosing the high density 
that matches this ratio under-estimates the other 



nebular line fluxes because their critical densites 



are much lower than that of [O III] AA5007,4959, 
which have critical density, n c = 6.5xl0 5 cm -3 . 

The detailed analysis of P84 showed this could 
be remedied by employing two separate gas com- 
ponents having different densities. In this sce- 
nario, the nebular lines overwhelm their auro- 
ral counterparts in a low density region where 
their ratios are determined solely by tempera- 
ture. Conversely, in a higher density gas the 
nebular line emissions are suppressed due to col- 
lisional de-excitation and, consequently, the au- 
roral line fluxes become a more significant source 
of cooling. Hence, a density stratification effect is 
mistaken for an indication of high temperature. 
Multi-component gases having a range of densi- 
ties are a common feature in models for LINERs 
and Seyfert galaxies (Halpern & Steiner 1984; 
Filippenko 1985; Kraemer et al. 1994, 1998). 
Since we seek the simplest possible model, our 
method was to find the fewest number of dis- 
crete cloud components that would reproduce 
the observed emission. Although assuming a 
small number of distinct isolated gas components 
might be an oversimplification, it provides insight 
into the physical conditions and structure of the 
emission line gas. An alternate approach, the 
locally optimally emitting cloud model, has also 
been shown to effectively model the narrow line 
region emission in AGNs ( Ferguson et al. 1997| ). 

Thus, as a first modification to the simple 
model, we assumed a gas comprised of two den- 
sity components. Our method was to fit the 
nebular lines that have very low critical den- 
sities ([S II] AA6716,6731, [N II] AA6548,6584, 
and [O II] A3727) with a low density gas, com- 
ponent 1, and the strongest auroral lines (ie, 
those of S + , + , and +2 ) with a higher den- 
sity gas, component 2. Critical densities of var- 
ious lines place constraints on the density for 
each component. If the [O II] A3727 line is a 
principal coolant for component 1, the gas den- 



sity must not be significantly greater than its 
critical density of n c = 4500 cm~ 3 ( DeRober- 
tis fc Osterbrock 1984| ). A lower limit on the 
density for component 1 is obtained from the 
[S II] AA6716,6731 doublet lines. The flux ra- 
tio [S II] A6716/ [S II] A6731 = 1.1 indicat es an 
electron density n e > 300 cm -3 (|Aller 1982j ), For 
the high density component, the [O III] A4363/ 
[O III] A5007 ratio implies 6.5 xlO 5 < n e < 
3.3xl0 7 cm -3 as determined by the critical den- 
sities of the two lines. Additionally, the strong 
[O I] A6300 emission indicates n e < 1.8x 10 6 cm -3 
QDeRobertis fc Osterbrock 1984| ). 

Since the spectrum of NGC 1052 is dominated 
by lines of low ionization species, lines emitted by 
S + , N + , O , + , and C + , a low ionization pa- 
rameter is required. Figure 1 of Ferland & Net- 
zer (1983), which shows the relationship between 
various line strengths and the ionization param- 
eter for a low density gas, indicates U < 10~ 3 ' 5 
for component 1. We find that a slightly higher 
U is needed for component 2 to sufficiently pop- 
ulate the states responsible for the auroral line 
emission. Model results using n# = 10 3 ' 5 cm -3 , 
Ui = 10~ 3 - 8 and n H = 10 6 2 cm" 3 , U 2 = 10 -31 
for components one and two, respectively, are 
presented in columns 3 and 4 of Table 4. The 
composite model in column 5 of Table 4 was cre- 
ated by assuming equal contribution from each 
component to the total observed H/3 flux. 



3.2. The Effects of Grains in the Emis- 
sion Line Gas 

Inspection of Table 4 shows that the overall 
intrinsic emission line spectrum of NGC 1052 is 
fit very well by our photoionization model, which 
has employed a minimum of assumptions. Partic- 
ulary, the strong auroral line fluxes of [O III] A4363, 
[O II] A2470, and [S II] AA4069,4076 are matched 
well, as are the [O III] A4363/[0 III] A5007 and 
[O II] A2470/ [O II] A3727 ratios, which are of- 
ten mistakenly used as temperature indicators. 
However, several of the prominent low ioniza- 
tion features are considerably underpredicted. 
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For example, the observed [N II] AA6548,6584, decrease in their strengths is expected with de- 

[S II] AA6716,6731, and [O I] AA6300,6363 line pletion. Indeed, comparison of the results of our 

fluxes are matched within a factor of two, but are two models demonstrates this. Specifically, the 

all somewhat underestimated. The C II] A2326 dusty model predicts a decrease in the strength of 

fit is even worse. the lines by a factor of ten. Although the [Ca II] 

All of these poorly fit lines emit predominately lin es lie outside of the spectral region covered by 

in the warm partially ionized zones, [S II] AA6716, °ur da ^ the near infrared spectral data of Diaz 

6731 and [N II] AA6548,6584 most strongly in et al. (1985) clearly show that their strengths are 

component 1 and [O I] AA6300, 6363 and C II] A2326 overpredicted by our dust free model (see their 

Figure 1). Thus, there is evidence that signifi- 
cant depletion exists in NGC 1052. 

Our depletion percentage corresponds to a gas 
column-to-reddening scaling of roughly Njy /E(B- 
V) = lxlO 22 cm" 2 mag" 1 flKnapp h Kerr 1974|) . 
The ionized column densities of our gas compo- 
nents are comparatively small (Nhii = 3.1 xlO 19 
cm -2 and Nhii = 3.6 xlO 20 cm~ 2 , for compo- 
nents 1 and 2, respectively). Thus, most of the 
reddening occurs in a neutral gas obscuring our 
line of sight to the emitting regions. If this ob- 
scuration is due to gas component 1 of our model, 
n# = 10 35 cm" 3 , a cloud thickness of only ~ g pc 
is required to produce the observed extinction. 

The polarimetry study by Barth et al.(1999) 
shows that [O I] A6300 is more polarized than the 
other "narrow" emission lines between ~4800— 
6800 A, with a polarization percentage more 
than twice as much as any other line. This 
indicates that [O I] A6300 may be transmitted 
through more dust than the other lines. This 
effect is explained naturally, as in our model, if 
[O I] A6300 is emitted predominately in an ex- 
tended dusty partially ionized zone (PIZ). If we 
are viewing the line emission primarily from the 
Ne = 8.51xl0- 5 ,Mg = 1.07xl0- 5 ,Si = 1.20xl0- 5 , ionized faces of the gas clouds, then the lines 
S = 1.29xl0~ 5 , and Fe = 1.02xl0~ 5 . emitted in the "backside" of the gas pass through 

Model results, using the same parameters as in a much larger dust column before escaping. Our 

our previous model are given in Table 5. We have calculations show that the PIZ of component 2, 

included model predictions of the [Ca II] AA7291, where [O I] A6300 is most strongly produced, has 

7324 lines for both the dusty and solar abundance a much larger column density than the PIZ of 

models in Tables 5 and 4, respectively. Ferland component 1. Since this is the region where the 

(1993) has shown that the weakness of these lines [N II] AA6548,6584 and [S II] AA6716,6731 lines 

in the narrow line regions of some AGNs is a are emitted, this explains why there is much less 

good indicator of the presence of grains - a sharp polarization of these lines than of the [O I] A6300 



primarily in the high density component 2. The 
characteristic emission in these regions is likely 
greatly modified by the presence of dust due to 
the heavy depletion of important coolants (par- 
ticularly magnesium and silicon) and an increase 
in gas temperature caused by the selective re- 
moval of photons near the Lyman limit by dust 
grains. To treat this, we have included grains 
typical of the Orion Nebula mixed with the emis- 
sion line gas ( Baldwin et al. 1991| ). We se- 
lected this type of grain structure since, of the 
available data for dust grains, physical condi- 
tions in an H II region most closely approximate 
those of NGC 1052. Constraints on the dust-to- 
gas ratio in NGC 1052 come from the strengths 
of the Si III] AA1882,1892 and Mg II A2800 
fluxes. The calculated emission of these lines 
becomes too weak when the full depletion de- 
duced for Orion is used in our models (Ferland 



19$t]: based on r esults of Baldwin et al. 1991 



Rubin et al. 1991 , and Osterbrock 1992 ). There- 
fore, we assumed half of the depletion and grain 
abundances, giving the following elemental abun- 
dances relative to hydrogen: He = 9.77xl0~ 2 , 
C = 3.31xl0~ 4 , N = 8.13xl0~ 5 , O = 5.50xl0" 4 , 
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line. 

3.3. Structure of the Emission Line Re- 
gion 

We now turn our attention toward the struc- 
ture of the emission line region as deduced from 
our model results. Since we have an estimate 
of the luminosity of ionizing radiation, we can 
calculate the radial distances of the different gas 
components from the putative central engine, Ri 
and R2, using U oc L^^^/R 2 . This gives 
Ri ~ 100 pc and R2 ~ 2 pc. The deduced lo- 
cation of component 1 lies outside the projected 
aperture by about a factor of 2, however, the 
uncertainty in the intrinsic luminosity could ac- 
count for this discrepancy. Furthermore, if the 
gas is spherically distributed, then much of it 
would be projected into the aperture. 

The broad widths of the emission lines (900 < 
FWHM < 1100 km s" 1 ) give clues to the gas dy- 
namics in NGC 1052. We find no correlation be- 
tween line width and critical density. Thus, since 
the high and low density components are at very 
different distances from the putative central mas- 
sive object, the velocities of the two components 
cannot be explained as simple circular Kepler ian 
orbits, for which v oc r -1 / 2 . 

As mentioned earlier, our assumption that the 
line emitting gas is comprised of distinct compo- 
nents having constant densities is a rough ap- 
proximation, but qualitatively the higher den- 
sity gas is seen to be closer to the center of the 
nucleus. To match the calculated H/3 luminos- 
ity with the observed luminosity, covering fac- 
tors of ei = 0.28 and ei = 0.25 are required for 
components 1 and 2, respectively, assuming that 
each component contributes equally to the total 
H/5 emission. Model predictions of ~Lh@ f° r each 
component are given in Table 5. Since the cov- 
ering factors for each component were calculated 
independently of each other, the attenuation of 
the inner component was not treated and hence 
ei represents a minimum covering factor for the 
outer component. Therefore, a minimum total 



covering factor for NGC 1052 is e*r = 0.53, which 
represents the probability that the line-of-sight to 
the central continuum is obscured by the emis- 
sion line gas. This is considerably higher than 
that determined for the narrow line regions of 
Seyfert 1 galaxies ( Netzer fc Laor 1993j ). Thus, 
our previous assumption that the central ioniz- 
ing continuum is occulted by a dusty gas is vin- 
dicated. This is consistent with the study by 
Barth et al. (1998) which concluded that many 
LINERs may have bright compact UV sources 
that are obscured by dust. 

3.4. The High Ionization Component 

Although the spectrum of NGC 1052 is dom- 
inated by low ionization features, there are sev- 
eral high ionization features present, namely, 
C IV A1549, [Ne IV] A2424, and [Ne V] A3426. 
As anticipated from our choice of low ionization 
parameters, Ui & U2, these lines are greatly 
underestimated. In particular, the predicted 
C IV A 1549 is ~4 times too weak and [Ne IV] A2424 
and [Ne V] A3426 are underpredicted by a factor 
of ~10. Imposing the required U on either of the 
two gas components would severely degrade our 
fit to the prominent low ionization features, giv- 
ing results characteristic of Seyfert galaxies. In- 
deed, it is this lower ionization that distinguishes 
LINERs from other active galaxies. Thus, there 
is an indication that a highly ionized gas compo- 
nent is present in NGC 1052, although its con- 
tribution is small. 

The critical density of [Ne IV] A2424 limits 
n e <3xl0 4 cm -3 . However, for a gas with this 
low density and having a high ionization param- 
eter, the [O III] AA5007,4959 lines become the 
dominant coolants and their fluxes are far over- 
predicted. We find that an emitting gas that is 
matter bounded can account for this. For illus- 
trative purposes, we have calculated the emission 
for a gas without dust grains having U3 = 10~ L5 
and iih = 10 4 2 cm -3 and a total column den- 
sity of N# = 10 21 cm -2 . This corresponds to a 
distance R3~3 pc of the highly ionized gas from 
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the central continuum source. A very low cov- 
ering factor of only ~0.05 is required to match 
the observed C IV A 1549 luminosity with this 
extra component, hence our assumption that all 
components receive unattenuated continuum ra- 
diation remains valid. This simple model for a 
highly ionized component matches the observed 
[Ne V] A3426 emission line very well and the 
[Ne IV] A2424 line within a factor of 2. Further- 
more, since the emissivity of C IV A1549 is very 
high under these conditions, there would be only 
a small contribution from this component and, 
therefore, no change in the low ionization line ra- 
tios predicted by a three-component model from 
those listed in Tables 4 & 5. 

4. Discussion 

We have shown that a simple photoionization 
model explains the full set of dereddened opti- 
cal and UV emission lines observed in the FOS 
data of NGC 1052 quite well. We note that the 
reddening correction we have applied to the UV 
continuum is highly uncertain and that this en- 
hanced luminosity is essential in reproducing the 
line luminosities. Our best fit model employs 
a dusty emission line gas with a high covering 
factor. This is consistent with the apparent ab- 
sorption of the soft X-ray flux observed in the 
ROSAT data, indicating that the line-of-sight to- 
ward the AGN is obscured. The flat intrinsic 
ionizing continuum (a = —1.2) that we have de- 
duced for NGC 1052 is in accord with the recent 
study by Ho (1999) which demonstrated that a 
sample of LINERs all have a higher X-ray/UV 
continuum ratio than found in typical high lumi- 
nosity AGNs. 

4.1. Extended Emission Line Region 

We examine if the deduced ionizing contin- 
uum flux is sufficient to power the extended emis- 
sion line region (EELR) emission described in 
F78. They found a diffuse emission line gas ex- 
tending out to about 20" from the nucleus hav- 



ing a total Balmer line flux similar to that ob- 
served in the nuclear emission region. Using 
their flux measurements and assuming that all 
of the emission line fluxes experience the extinc- 
tion that we have deduced from our analysis gives 
L# ( g~2.6x 10 40 ergs s _1 for the total H/3 luminos- 
ity in the extended plus compact emission com- 
ponents. A simple photon-counting argument 
( Netzer fc Laor 1993|) using our deduced contin- 
uum shape and Lyman edge flux demonstrates 
that a covering factor of ~2 is required to pro- 
duce the H/3 luminosity. 

A covering factor greater than unity indicates 
that if all of the line emission in NGC 1052 (ex- 
tended + compact) is powered by photoioniza- 
tion from the putative non-thermal continuum, 
then our view of this continuum must be even 
more obscured than we have previously deduced. 
Furthermore, the large covering factor of the in- 
ner gas regions calculated in our model suggests 
that the shape of the ionizing continuum inci- 
dent upon the extended region would be greatly 
modified. However, the attenuation of UV con- 
tinuum photons by the emission line gas in the 
compact region would result in a much harder 
SED incident upon the extended gas and a much 
higher actual covering factor for the EELR is re- 
quired. An alternate explanation is that the cen- 
tral source emission is beamed in NGC 1052 such 
that the continuum flux "seen" by the emitting 
gas in the extended regions is enhanced. An- 
other possibility, is that some process other than 
photoionization by the central source drives the 
emission in the EELR, either stellar photoioniza- 
tion or shock excitation. A detailed analysis of 
the diffuse emission regions is needed to fully test 
this. 

4.2. Effects of the Choice of Extinction 
Curve 

As noted in section 3, one discrepancy in our 
model is the underestimate of the C II] A2326 flux 
by a factor of ~2. The C II] A2326 emission is 
strong in a gas with low ionization parameter and 
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high density. However, calculations show that 
the Mg II A2800 flux is far overpredicted using 
these model conditions. Alternately, C II] A2326 
is produced predominantly in the PIZ of a pho- 
toionized gas, suggesting a stronger X-ray con- 
tinuum would strengthen the feature. However, 
the [O I] AA6300,6363 line fluxes are more sensi- 
tive to the high energy flux and are greatly over- 
estimated if the X-ray continuum is enhanced. 

A poor choice of the extinction curve may 
account for this discrepancy. Several studies 
have shown that the UV dust extinction in ex- 
tragalactic objects is not represented well by 



the Galactic Savage & Mathis curve ( Fitzpatrick 
19^, |Calzetti et al. 1994| ). Often a weaker 



2200 A absorption feature and steeper rise in the 
far-UV are found. As an ad hoc demonstration, 
we show how the UV line strengths are affected 
by removing the 2200 A bump from the extinc- 
tion curve in Table 6. The [C II] A2326 and 
[Ne IV] A2424 line fluxes are both significantly 
reduced, illustrating that if an extinction law in- 
corporating a strong 2200 A feature were erro- 
neously applied, it could lead to a substantial 
overestimate of the intrinsic strengths of these 
lines. 

5. Conclusions 

We have analyzed archival HST/FOS data of 
the prototypical LINER galaxy NGC 1052. The 
aperture sampled the inner ~ 82 pcx82 pc of the 
nuclear region, thereby excluding the extended 
emission region. One important result is that the 
deduced Ha/H/3 ratio provides evidence for con- 
siderable extinction (E(B-V) =0.42) of the emis- 
sion line fluxes in NGC 1052. Dust appears to 
be concentrated in or near the nucleus and may 
also be attenuating the nonthermal continuum 
of a compact source, which supports the Barth 
et al. (1998) conclusion that the nuclei of many 
LINERs are obscured by dust. 

From our modeling results, we find that the 
nuclear emission observed in the inner region of 



NGC 1052 is explained well by photoionization 
from a simple power law continuum having a 
flat spectral index, a=— 1.2. Moreover, the ion- 
izing continuum is consistent with our deduced 
reddening and observed flux distribution. This 
is also consistent with the study by Ho (1999) 
which showed that a sample of LINERs all have 
a flat spectral energy distribution in the ionizing 
continuum. 

Most of the dereddened line fluxes are matched 
by a model having two dusty gas components 
of different densities, photoionized by a single 
power law continuum. The appearance of the 
C IV A1549, [Ne IV] A2424, and [Ne V] A3426 
lines in the spectrum indicates a third gas com- 
ponent may be present, characterized by a high 
ionization parameter. The key exception to our 
good model fit is the prominent C II] A2326 line, 
which is underpredicted by a factor of about 
two. This flux is heavily affected by the strength 
of the 2200 A absorption feature. Hence, the 
C II] A2326 line would be overestimated if the 
2200 A bump of the Savage & Mathis extinc- 
tion curve is not applicable to NGC 1052. This 
would also help explain the underestimate of the 
[Ne IV] A2424 flux. 

The large widths of the emission lines and lack 
of correlation between width and critical density 
suggests non-gravitational motion of the emission 
line gas. 

We conclude that a pure central source pho- 
toionization model with the simplest non-thermal 
continuum (a simple power law) reproduce the 
emission line fluxes in the inner region of NGC 1052 
quite well. Other processes such as shocks or 
photoionization by stars are not required to pro- 
duce the observed emission. However, we cannot 
rule out the possibility that these mechanisms 
contribute to the extended nebular emission. 

This research has made use of the SIMBAD 
database operated at CDS, Strasbourg, France 
and of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion 
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Fig. 1. — FOS spectra of the nucleus of 
NGC 1052 from 1250 to 6800 A, corrected for 
Galactic reddening. The plotted template spec- 
tra are scaled to match the absorption features 
of NGC 1052 and offset in flux for clarity. 

Fig. 2.— (2a) Deblending procedure for 
Ha+[N II]. The solid line is the Ha+[N II] line 
blend with the host galaxy template removed. 
The [N II] AA6548,6584 emission line template 
(dot-dashed line) was created from the H/3 pro- 
file and subtracted to give the Ha (dashed line) 
emission line. (2B) Comparison of the shape of 
the deblended Ha profile to the H/3 profile scaled 
to match in flux. 

Fig. 3. — The multiwavelength spectrum of the 
nucleus of NGC 1052. The radio data (+) are the 
measurements of the compact emission (< 1") 
by Wrobel et al. 1984. The IR data points (*) 
are within < 2" (Becklin et al. 1982). The ob- 
served and dereddened UV continuum data at 
1300 A are from this study ~0".86 (o). ASCA 
and ROSAT data (A) are from Guianazzi & An- 
tonelli (1999) with a resolution of 5'. It is unclear 
how much of this X-ray emission is compact. The 
ionizing continuum used in our photoionization 
models is shown as the solid line. The dashed 
line is the attenuated continuum, assuming the 
central engine is viewed through the high density 
gas, component 2. 

Fig. 4.— WFPC2 F658N image of the inner 
region of NGC 1052. The size of the FOS 0."86 
square aperture is plotted over the compact emis- 
sion core for reference. North is up in this image. 



16 



Table 1 

FOS Observations of NGC 1052 



Data Set 


Grating 


Resolution 

(A) 


Exposure 
(s) 


Template Galaxy 


Y3K80105T 


G160L a 


-6.4 b 


1390 




Y3K80106T 


G160L 




2430 




Y3K80107T 


G160L 




2430 




Y3K80108T 


G270H 


~2.1 c 


2320 


NGC 3608 


Y3K80109T 


G400H 


~3.1 c 


1200 


NGC 3608 


Y3K8010AT 


G570H 


~4.4 C 


890 


NGC 5845 



a G160L spectra coadded, weighted by exposure time 

b A/AA~250 

C A/AA~1300 
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Table 2 

Measured and Dereddened Emission Line Ratios' 1 





E(B-V)= 


=0.02 


E(B-V> 


=0.44 


C TV A154Q 


0.14 


(± 


04) 


0.75 


V 


64) 


Hp TT A1 640 

-L-L /\ lU^U 


0.16 


(± 


05) 


0.80 


V 


66) 


ro TTTl A1663 


0.07 


(± 


02) 

\J • UZv \ 


0.35 


V 


28) 


N TTTl A1 75D 


0.04 


(± 


n m l 


0.20 


(± 

V 


n 1 6) 

\J. -L\J I 


Q\ TT 4. [Wp TTTl A1 808-1 820 


0.10 


(± 


03) 








Si TTTl AA1882 1 8Q2 b 


0.14 


(± 


n 05) 

yj.vjo } 


0.80 


(± 


n 71 ) 


C TTTl A1909 b 


0.35 


(± 


08) 


2.00 


(± 


1 73) 

L. 1 U 1 


N II] A2140 


0.06 


(± 


0.01) 








C II] A2326 c 


0.58 


(± 


0.08) 


3.80 


(± 

V 


3.51) 


[Ne IVl A2424 

J. i C -L V 


0.02 


(± 


0.01) 


0.10 


(± 


0.09) 


O III A2470 

V/ -L-L 1 \J 


0.14 


(± 


0.03) 


0.60 


(± 


0.43) 


Fe II A2558-2583 


0.03 


(± 


0.01) 








Fe II A2601-2640 


0.14 


(± 


0.03) 








Fp TT A2660-2684 

_L C -L-L /\ Zj \J\J\J Zj \J Ot: 


0.04 


(± 


01 ) 








Fp TT A271 1-2765 

J. C- -L-L /\ 1 X -L ( *J 


0.15 


(± 


03) 








Fe II A2823-2844 


0.04 


(± 


0.01) 








Mo- TT A98Dn 

IV-Lg _L_L /\ZvOL/VJ 






•fl 46 d 






1 2 


Hp TT A3204 

1 1 V -L-L /\ kJ j-j\J^ 






••0 05 






••0 08 


[Ne VI A3346 






0.015 






::0.03 


Ne V A3426 


0.07 


(± 

V 


0.02) 


0.11 


V 


0.04) 


[0 II] A3727 


2.77 


(± 

V 


0.31) 


3.98 


V 


0.81) 


[Ne TTTl A 3869 


0.60 


(± 


08) 


0.82 


(± 


1 6) 


Np TTTl AS967 b 


0.20 


(± 


n ns) 


0.27 


(± 


n 05) 


He A3970 b 

1 1 V / \KJ \J 1 V / 


0.14 


(± 


0.03) 


0.19 


(± 

V 


0.05) 


FS III AA4069 4076 


0.68 


V 


0.10) 


0.88 


(± 


0.17) 


H<5 A4102 


0.24 


(± 


0.03) 


0.31 


(± 


0.05) 


H^ A4340 b 

-L-L J AltjT:U 


0.44 


(± 


08) 


0.53 


(± 


11) 


[0 III] A4363 b 


0.27 


(± 


0.05) 


0.32 


(± 


0.06) 


He II A4686 






::0.18 






::0.19 


H/3 A4861 






1.00° 






1.00 


[0 III] A4959 


0.91 


(± 


0.13) 


0.88 


(± 


0.13) 


[0 III] A5007 


2.85 


(± 


0.32) 


2.72 


(± 


0.31) 


[N I] A5200 


0.23 


(± 


0.03) 


0.21 


(± 


0.03) 


[N II] A5755 






::0.09 






::0.07 


He I A5876 






::0.14 






::0.10 


[0 I] A6300 


3.32 


(± 


0.37) 


2.22 


(± 


0.49) 


[0 I] A6364 


1.13 


(± 


0.13) 


0.74 


(± 


0.17) 


[N II] A6548 b 


1.73 


(± 


0.26) 


1.10 


(± 


0.29) 


Ha A6563 b 


4.53 


(± 


0.68) 


2.86 


(± 


0.77) 


[N II] A6583 b 


is 5 - 14 


(± 


0.77) 


3.24 


(± 


0.87) 



Table 2 — Continued 





E(B-V)=0.02 


E(B-V)=0.44 


He I A6678 


::0.07 


::0.04 


[S II] A6716 b 


2.56 (± 0.38) 


1.57 (± 0.43) 


[S II] A6731 b 


2.35 (± 0.35) 


1.45 (± 0.40) 



a Relative to H/3 

b Line deblending procedure required for flux 
measurement 

c Contains [O III] A2321 and Si II] A2335 
emission 

d Mg II A2800 is heavily absorbed; a fit 
to the doublet with Gaussian features gives 
Mg II A2800/H/3 = 0.80 

e F Hf} = 2.65(±0.13)xl0" 14 ergs s" 1 cm" 2 

f F^ = l.l(±0.7)xl0- 13 ergs s" 1 cm" 2 
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Table 3 

Intrinsic and Dereddened Balmer Line Ratios 3, 





Intrinsic b 


Dereddened 


Ha A6563 


2.86 


2.86 (± 0.77) 


H 7 A4340 


0.47 


0.53 (± 0.11) 


H<5 A4102 


0.26 


0.31 (± 0.05) 


He A3970 


0.16 


0.19 (± 0.05) 



a Relative to H/3 

b From Osterbrock 1989 
c E(B-V) T =E(B-V) G +E(B-V)/=0.44 
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Table 4 

Line Ratios From Model without DusT a 





Observed 13 


Comp l c 


Comp 2 d 


Composite 6 


C, TV A 1 ^40. 

w IV AlO^iy 


(17^4(1 RA\ 




n 98 


n 14 

U.14: 


Hp tt m fi/in 

lie 1J. A1U4:U 


n on (-L. o p.a\ 
u.ou im u.uu I 


1 nQ 


1 98 


1 1Q 
i.iy 


[O TTTl Mfifi^ 

U lllj A1UUO 


n (-\-0 oq\ 
u.oo tin u.iu I 




n 4Q 




N TTTl A1 75n 


n on C+ f) Ifil 




90 
u. z>u 


n in 

u. ±u 


TTTl \1 889 1 8Q9 


n an (4- n 71 
U.oU ^zn u. / 1J 


n m 

U.Ul 


1 1 9 
1 . 1Z 


n ^7 

U.O i 


P TTTl \1 QDQ 


o nn ('-i- i 7*^ 

z.uu ^zn l. / oj 


n 9i 

U.Zl 


^n 
o.ou 


1 8fi 
l.ou 


n TTl 

±±J AZOZU 


o.ou o.oij 


i nn 

l.UU 


z.ou 


1 78 
1 . / o 


[M„ TVl \9/19/| 
[IN e 1 V J AZ41Z41 


n i n /'-i- n na^i 
u.iu ^zsz u.uyj 




n ni 

U.Ul 


n ni 

U.Ul 


rn TTl A247D 


n fif) f+ n 43^1 


30 

u.ou 


1.11 


71 

U. 1 -L 


1\/T<r TT \9ann 
ivig ii azouu 


••1 9 
..l.Z 


1 77 
1. 1 1 


a ^n 


^! 1 A 
0. 141 


TTo TT X^n/L 

lie 11 AOZU^fc 


• •n ns 

. .U.UO 


n n7 
u.u / 


n ns 

U.Uo 


n ns 

U.Uo 


HNFp VI X^df, 

llNe VI AOO^iU 


••n 

. .U.UO 






<rV\ ni 


I"Np VI A^i49fi 

I1NC VI AO^l^U 


nil f _i_ n fizL^l 
u.ii izn u.u^ii 






^n ni 


TO TTl X3797 

[vj 11J AO / Z / 


q no a nsi'i 
o.jo ^zn u.ol J 


7 n7 

i .U / 


n na 
u.uy 


^ ^8 
O.Oo 


riVo TTTl X^SfiQ 
[inp iiij Aoooy 


n Q9 (-1. c\ -\a\ 


n 89 

U.OZ 


1 Q9 
i.yz 


1 ^7 
1.0 / 


TlVTo TTTl \":!QR7 
[IN e lllj AoaD i 


n 97 / -i- n r\K\ 

U.Z/ ^It U.UOJ 


n 9^ 
u.zo 


n Rn 

U.DU 


n /i^ 

U.4o 


TT^ A^Q7D 


n in (4-n r\K\ 
u.iy tin u.uot 


n 1 7 

U. 1 / 


n 1fi 

U.IU 


n 17 

U. 1 1 


TCI TTl \ \4flRQ zL(T7fi 
[o 11J AA4UOy ,4LU / 


U.oo ^zn U. 1 1 J 


n zLQ 
u.^ty 


n Q9 
u.yz 


n 7i 

U. / 1 


TTA \zL1fl9 
no A^iuz 


n 01 /_|_ n nc^ 


n 97 

u.z / 


n 97 

U.Z ( 


n 97 

U.Z ( 


TT^v XzL^zLfl 

A4L041:U 


n ^ (A- n 1 1 
u.oo \pi u. 1 1 ) 


n 47 


n 47 


n 47 


rn ttti Xzt^fi^ 

[W lllj A^IOUO 


n Q9 (-L. n nfi"l 
u.oz u.uuj 




n 

u.oy 


n ^n 

u.ou 


TTo TT \/1fi8f; 
116 11 A4uo0 


• •n i q 
. .u. iy 


n 1 9 
u. 1Z 


n i A 

U. 14: 


U. lO 


TT/9 


1 nn 

l.UU 


i nn 

l.UU 


1 nn 

l.UU 


1 nn 

l.UU 


TO TTTl A4Q5Q 


n qo fj.niQ) 

U.OO I _l_ u.lO 1 


n 34 


1 79 


1 03 


O TTTl A 5007 


o 79 (a- n o9\ 


n Q7 


4 98 

t: . i?0 


9 98 

. i?0 


[N I] A5200 


0.21 (± 0.03) 


0.25 


0.01 


0.13 


[N II] A5755 


::0.07 


0.05 


0.17 


0.11 


He I A5876 


::0.10 


0.17 


0.15 


0.16 


[0 I] A6300 


2.22 (± 0.49) 


1.23 


2.19 


1.71 


[0 I] A6364 


0.74 (± 0.17) 


0.39 


0.70 


0.55 


[N II] A6548 


1.10 (± 0.29) 


1.36 


0.29 


0.83 


Ha A6563 


2.86 (± 0.77) 


2.92 


2.96 


2.94 


[N II] A6583 


3.24 (± 0.87) 


4.03 


0.85 


2.44 


He I A6678 


::0.04 


0.05 


0.03 


0.04 


[S II] A6716 


1.57 (± 0.43) 


1.47 


0.06 


0.77 


[S II]A6731 


1.45 (± 0.40) 


2.03 


0.14 


1.09 


[Ca IIJA7291 




0.61 


0.06 


0.34 


[Ca IIJA7324 




0.40 


0.04 


0.22 



Relative to H/3 21 

b Corrected for Galactic and internal reddening, E(B-V)=0.44 
C U=10~ 3 - 8 , N//=10 3 - 5 



Table 5 

Best-Fit Model:Includes Dust Grains a 



Observed 13 Comp l c Comp 2 d Composite 



C IV A1549 





75 


(± 


0.64) 







38 





19 


He II A1640 





80 


(± 


0.66) 


0.84 


1 


28 


1 


06 


[0 III] A1663 





35 


(± 


0.28) 


0.01 





53 





27 


N III] A1750 





20 


(± 


0.16) 







25 





13 


Si III] AA1882,1892 





80 


(± 


0.71) 


0.02 





54 





28 


C III] A1909 


2 


00 


(± 


1.73) 


0.35 


4 


53 


2 


44 


C II] A2326S 


3 


80 


(± 


3.51) 


1.45 


2 


86 


2 


16 


[Ne IVl A2424 





10 


(± 


0.09) 







01 





01 


[0 II] A2470 





60 


(± 


0.43) 


0.36 


1 


02 





69 


Mg II A2800 








::1.2 


0.89 


1 


16 


1 


03 


He II A3204 








::0.08 


0.05 





08 





07 


[Ne V] A3346 








::0.03 








<0 


01 


[Ne V] A3426 





11 


(± 
v- 1 - 


0.04) 








<0 


01 


[0 II] A3727 


3 


98 


(± 


0.81) 


7.35 





09 


3 


72 


[Ne III] A3869 





82 


(± 


0.16) 


0.78 


1 


68 


1 


23 


[Ne III] A3967 





27 


(± 


0.05) 


0.23 





51 





37 


He A3970 





19 


(± 


0.05) 


0.17 





16 





16 


[S II] AA4069,4076 





88 


(± 


0.17) 


0.55 


1 


06 





81 


H<5 A4102 





31 


V 


0.05) 


0.26 





27 





27 


H7 A4340 





53 


(± 


0.11) 


0.47 





47 





47 


[0 III] A4363 





32 


(± 


0.06) 


0.01 





58 





30 


He II A4686 








::0.19 


0.09 





14 





12 


H/3 h 

r 








1.00 


1.00 


1 


00 


1 


00 


[0 III] A4959 





88 


(± 


0.13) 


0.33 


1 


50 





92 


[0 III] A5007 


2 


72 


(± 


0.31) 


0.96 


4 


34 


2 


45 


[N I] A5200 





21 


(± 


0.03) 


0.40 





02 





21 


[N II] A5755 








::0.07 


0.06 





19 





13 


He I A5876 








::0.10 


0.16 





14 





15 


[0 I] A6300 


2 


22 


(± 


0.49) 


1.40 


2 


49 


1 


95 


[0 I] A6364 





74 


(± 


0.17) 


0.45 





79 





62 


[N II] A6548 


1 


10 


(± 


0.29) 


1.44 





36 





90 


Ha A6563 


2 


86 


(± 


0.77) 


2.93 


2 


92 


2 


93 


[N II] A6583 


3 


24 


(± 


0.87) 


4.26 


1 


05 


2 


66 


He I A6678 








::0.04 


0.04 





03 





04 


[S II] A6716 


1 


57 


(± 


0.43) 


1.72 





10 





91 


[S II]A6731 


1 


45 


(± 


0.40) 


2.20 





23 


1 


22 


[Ca II]A7291 










0.07 





01 





04 


[Ca II]A7324 










0.04 





01 





03 



a Relative to H/3 22 

b Corrected for Galactic and internal reddening, E(B-V)=0.44 
C T7=10~ 3 - 8 , N H =10 3 - 5 



Table 6 

The Effect of the 2200 A Feature in the Extinction Curve 



S & M Curve a Modified S & M Curvef Best-Fit Model 



C IV A1549 


0.75 


(± 


0.64) 


0.75 


0.19 


He II A1640 


0.80 


(± 


0.66) 


0.80 


1.06 


[0 III] A1663 


0.35 


(± 


0.28) 


0.35 


0.27 


N III] A1750 


0.20 


(± 


0.16) 


0.20 


0.13 


Si III] AA1882,1892 


0.80 


(± 


0.71) 


0.70 


0.28 


C III] A1909 


2.00 


(± 


1.73) 


1.69 


2.44 


C II] A2326 c 


3.80 


(± 


3.51) 


2.01 


2.16 


[Ne IV] A2424 


0.10 


(± 


0.09) 


0.07 


0.01 


[0 II] A2470 


0.60 


(± 


0.43) 


0.39 


0.69 



a E(B-V) T = 0.44 

b Savage & Mathis extinction curve with 2200 A feature removed, E(B-V)t = 0.44 
c Contains [O III] A2321 and Si II] A2335 emission 



23 



F„ (10~ 15 erg s" 1 cm " 2 A" 1 ) 




F„ ( 1 O" 15 erg s" 1 cm " s A" 1 ) 




F„ (10~ 15 erg s" 1 cm " 2 A" 1 ) 




Normalized Flux (ergs s 1 cm 2 ) 




](!<>•( F ) (ergs s 1 cm 2 ) 




This figure "Fig4.jpg" is available in "jpg" format from: 
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